An improved effective-one-body Hamiltonian for spinning black-hole binaries 
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Building on a recent paper in which we computed the canonical Hamiltonian of a spinning test 
particle in curved spacetime, at linear order in the particle's spin, we work out an improved effective- 
one-body (EOB) Hamiltonian for spinning black-hole binaries. As in previous descriptions, we 
endow the effective particle not only with a mass fi, but also with a spin S,. Thus, the effective 
particle interacts with the effective Kerr background (having spin SKcrr) through a geodesic- type 
interaction and an additional spin-dependent interaction proportional to S,. When expanded in 
post-Newtonian (PN) orders, the EOB Hamiltonian reproduces the leading order spin-spin coupling 
and the spin-orbit coupling through 2.5PN order, for any mass-ratio. Also, it reproduces all spin- 
orbit couplings in the test-particle limit. Similarly to the test-particle limit case, when we restrict the 
EOB dynamics to spins aligned or antialigned with the orbital angular momentum, for which circular 
orbits exist, the EOB dynamics has several interesting features, such as the existence of an innermost 
stable circular orbit, a photon circular orbit, and a maximum in the orbital frequency during the 
plunge subsequent to the inspiral. These properties are crucial for reproducing the dynamics and 
gravitational-wave emission of spinning black-hole binaries, as calculated in numerical relativity 
simulations. 

PACS numbers: 04.25.D-, 04.25. dg, 04.25.Nx, 04.30.-w 



I. INTRODUCTION 

Coalescing black-hole binaries are among the most 
promising sources for the current and future laser- 
interferometer gravitational-wave detectors, such as the 
ground-based detectors LIGO and Virgo [l|, [Ij and the 
space-based detector LISA 

The search for gravitational waves from coalescing bi- 
naries and the extraction of the binary's physical param- 
eters are based on the matched filtering technique, which 
requires accurate knowledge of the waveform of the in- 
coming signal. Because black holes in general relativity 
are uniquely defined by their masses and spins, the wave- 
forms for black-hole binaries on a quasi-circular orbits 
depend on eight parameters, namely the masses mi and 
TO2 and the spin vectors 5'i and 82- Due to the large pa- 
rameter space, eventually tens of thousands of waveform 
templates may be needed to extract the gravitational- 
wave signal from the noise, an impossible demand for 
numerical-relativity alone. Fortunately, recent work at 
the interface between analytical and numerical relativity 
has demonstrated the possibility of modeling analytically 
the dynamics and the gravitational- wave emission of coa- 
lescing non-spinning black holes, thus providing data an- 
alysts with analytical template families [3-0] to be used 
for the searches (see also Ref. which considers the 
cases of extreme mass-ratio inspirals). The next impor- 
tant step is to extend those studies to spinning precessing 
black holes. 

So far, the analytical modeling of the inspiral, plunge ^, 



merger ^, and ringdown has been obtained within either 
the effective-one-body (EOB) formalism 0, S 0, H^ll or 
in Taylor-expanded PN models [l^, both calibrated to 
numerical-relativity simulations, or in phcnomenological 
approaches [l^ where the numerical-relativity wave- 
forms are fitted to templates which resemble the PN ex- 
pansion, but in which the coefficients predicted by PN 
theory are replaced by many arbitrary coefficients. Con- 
sidering the success of the EOB formalism in under- 
standing the physics of the coalescence of non-spinning 
black holes and modeling their gravitational-wave emis- 
sion with a small number of adjustable parameters, in 
this paper we will use that technique, adapting it to the 
case of spinning black- hole binaries. 

The first EOB Hamiltonian which included spin effects 
was computed in Ref. [l^ . In Ref. [2^ , the authors used 
the non-spinning EOB Hamiltonian augmented with PN 
spin terms to carry out the first exploratory study of the 
dynamics and gravitational radiation of spinning black- 
hole binaries during inspiral, merger and ringdown. More 
recently, Ref. [2l| extended the model of Ref. [19j to in- 
clude the next-to-leading-order spin-orbit couplings. The 
EOB formalism developed in Refs. [l^,[lH highlights sev- 
eral features of the spinning two-body dynamics and was 
recently compared to numerical-relativity simulations of 
spinning non-precessing black holes in Ref. [l^]. In this 
paper we build on Refs. [l^ [2l| and also on Ref. [l^l, in 
which we (in collaboration with Etienne Racine) derived 
the canonical Hamiltonian for a spinning test-particle in 
curved spacetime, at linear order in the particle's spin, 



the motion is driven mostly by the conservative dynamics. 
^ We refer to plunge as the dynamical phase starting soon after the ^ We refer to merger as the dynamical phase in which the two- 
two-body system passes the last stable orbit. During the plunge body system is described by a single black hole. 
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and work out an improved EOB Hamiltonian for spinning 
black-hole binaries. In particular, our EOB Hamiltonian 
reproduces the leading order spin-spin coupling and the 
spin-orbit coupling through 2.5PN order, for any mass- 
ratio. Also, it resums all the test-particle limit spin-orbit 
terms. Moreover, when restricted to the case of spins 
aligned or antialigned with the orbital angular momen- 
tum, it presents several important features, such as the 
existence of an innermost stable circular orbit, a photon 
circular orbit, and a maximum in the orbital frequency 
during the plunge subsequent to the inspiral. All of these 
features are crucial for reproducing the dynamics and 
gravitational- wave emission of spinning coalescing black 
holes, as calculated in numerical relativity simulations. 

This paper is organized as follows. After presenting 
our notation (Sec. HIl, in Sec. |lll] we build on Ref. [2J| 
and derive the Hamiltonian for a spinning test particle in 
axisymmetric stationary spacctimes. In Sec. IIVI we spe- 
cialize the axisymmetric stationary spacetime to the Kerr 
spacetime in Boyer-Lindquist coordinates. In Sec. IVl we 
work out the EOB Hamiltonian of two spinning precess- 
ing black holes. In Sec. I VII we restrict the dynamics to 
spins aligned or antialigned with the orbital angular mo- 
mentum and determine several properties of the circular- 
orbit dynamics. Section IVIII summarizes our main con- 
clusions. More details on how the spin-spin sector of the 
EOB Hamiltonian is constructed are eventually given in 
Appendix O 



II. NOTATION 

Throughout this paper, we use the signature 
(—,+,+,+) for the metric. Spacetime tensor indices 
(ranging from to 3) are denoted with Greek letters, 
while spatial tensor indices (ranging from 1 to 3) are 
denoted with lowercase Latin letters. Unless stated oth- 
erwise, we use geometric units (G = c = 1), although we 
restore the factors of c when expanding in PN orders. 

We define a tetrad field as a set consisting of a timelike 
future-oriented vector and three spacelike vectors ij 
(/ = 1, 3) — collectively denoted as (A = 0, 3) 



satisfying 



9 til/ 



VAB ■ 



(2.1) 



where titt = —1, Vti ~ 0, r//j = 6ij {6ij being the 
Kronecker symbol). 

Internal tetrad indices denoted with the uppercase 
Latin letters A, B, C and D always run from to 3, while 
internal tetrad indices with the uppercase Latin letters 
/, J, K and L, associated with the spacelike tetrad vec- 
tors, run from 1 to 3 only. The timelike tetrad index is 
denoted by T. 

Tetrad indices are raised and lowered with the metric 
rjAB [e.g., = r]AB (e'^)^]. We denote the projections 
of a vector V onto the tetrad with = ef: , and 
similarly for tensors of higher rank. Partial derivatives 



will be denoted with a comma or with 9, and covariant 
derivatives with a semicolon. 



III. HAMILTONIAN FOR A SPINNING 
TEST-PARTICLE IN AXISYMMETRIC 
STATIONARY SPACETIMES 

Following Ref. [23 |. we write a generic axisymmetric 
stationary metric in quasi-isotropic coordinates as 

ds^ = -e^^df + i?2 sin^ OB^e-^'' {d(f> - cudtf 

+6^^ (di?2 + i?^rf6|2) , (3.1) 

where v, fj,, B and uj are functions of the coordinates R 
and 9. Introducing the cartesian quasi-isotropic coordi- 
nates 



X = R sin 9 cos cj) , 
Y = R sin 9 sin <j> , 
Z = Rcos0, 

we can write Eq. (|3.1[) as 

ds2 ^ [^2^2 (^2 ^ - e"''] df 

+2B^ e-^" uj{YdX -X dY) dt 

-2-!^ dX dY 

X^ + Y^ 

g2M^2^52g-2.y2 



(3.2a) 
(3.2b) 
(3.2c) 



X^ + Y^ 
52g-2,. X^ + e^i^ y2 

X^ + Y^ 



dX' 
dY^ + e^^ dZ^ 



(3.3) 



It is straightforward to see that in the flat-spacetime 
limit {uj = ly = fj, = 0, B = 1) Eq. (|3.3p reduces to 
the Minkowski metric. 

Reference (23j computed the Hamiltonian of a spinning 
test-particle in curved spacetime at linear order in the 
particle's spin, and showed that it can be written as 



H — Hj<ss + Hi 



(3.4) 



where H^s is the Hamiltonian for a non-spinning test 
particle of mass m, given by 



H^s = P'P^ + aJm■2 + P^P^^ (3.5) 



with 



1 



ti 



9_ 

5* 



(3.6) 
(3.7) 
(3.8) 
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and 



Hs = -\ Fl" + Fl' 



a7^J Pi F^ 



Sk 



where the coefficients can be expressed in terms of a 
reference tetrad field ca as 



7K 



E 



LOT 



with 



Pt 
LOT 



P, 



P'Pi-a Jm'^ + ^3 P, Pj 



fL '^T — M T 



e: 

Lo,, e 



P 



(3.9) 
(3.10) 



(3.11) 
(3.12) 

(3.13) 

(3.14) 
(3.15) 



Reference [2j| also showed that in order to obtain a 
Hamiltonian giving the usual leading-order spin-orbit 



J 



coupling without gauge effects (or, equivalent ly, Hg — 
in flat spacetime) , the reference tetrad field must become 
cartesian in the flat-spacetime limit. We find that the 
following choice for the reference tetrad 



^1 



(J V H 



(3.16a) 



XY 



B{X^+Y^) 

(3.16b) 

e'')XY X^ + Be-t'Y^ 



B{X^ + Y^) ^ 
(Be-^ 



B (X2 + y2) 



■ A" 

B{X-^ + Y^) 

(3.16c) 
(3.16d) 



indeed reduces to the cartesian tetrad = 1, ef = 5f 
in the flat-spacetime limit. 

We can then use the tetrad defined by Eqs. (|3.16ap - 
(j3.16dp to calculate the coefficients F^ in Eq. p.9p . and 
obtain 



Hq — H[ 



so + ^^ss 



(3.17) 



with 



SO 



(e^+'^ ~B){P- ^R) 



B' 



R{S-0 YiR{P-yR)[VQ + ^)- A^cose {P-N)e-VQ {^R (P ■VR) + {p.^se - i^cose){P ■N)e) 
■e^+'iP-iR) (2v/0 + l) [iyRRiS-V)-Vcos9{S-N)e]B-BRe''+^P-^R) (^Q + l) R{S-V)j , 

(3.18) 



ifss =oj S 



e-^f'-'^LOR 



6^'+'' {p.VR){P-^R){S-i)B + e2(^+'^) (P • I R)^{S ■ V) 



2B iVQ + i) VQRe ^ 

e^" (i + Vq) ^/QR\S ■V)^'^ B^ + {P ■N)R {P-VR){S-N)-{P-N)R{S-V) 

(P • N){P ■iR)R{S-$,)B- e2(^+'')(P -iPfiS ■ N) 



e i^cos e 



2B + VO^' ^ 

{S■N){P■VR)^-{P■N)R{S■V){P■VR)-e'^^' (i + ^/q)^/QR'^{S-N)^^] B^\ , (3.19) 



Q = l+ Y^P^ Pj=l + (P • Nf' 



-2M(p.Vi?)2 e2'^(P-^i?)2 



i?2 ^2 



52 i?2 ^2 



(3.20) 



where we denote 



and 



ez X N ■ 



-Y ex +XeY 



(3.21) 
(3.22) 
(3.23) 



_ df{R,cos 
JR = ■ 



dR 
df{R, cost 
d{cos 9) 



(3.25) 
(3.26) 
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Here, the generic function / can stand for B, w, /i or v. 
Note that because uj is proportional to gt^ [see Eq. p.ip ] 
and thus to the spin of the spacetime, ifss (which is pro- 
portional to Lo and its derivatives) gives the leading-order 
coupling between the particle's spin and the spin of the 
background spacetime (together with other higher order 
terms). Also, because P • ^ i? = in spherical coordi- 
nates, i?so is the part of the Hamiltonian which gives the 
leading-order spin orbit coupling (again, together with 
other higher order terms). Moreover, note that iJs = 
in a flat spacetime, thus confirming the absence of gauge 
effects in the leading order spin-orbit coupling. 

As a consistency test, we specialize to the case of a 
spherically symmetric spacetime in quasi-isotropic coor- 
dinates, which was considered in Rcf. [2^ (see Sec. V 
A therein). Because the metric for such a spacetime is 
given by 

ds^ = -f(R) + h{R) {dX^ + dY^ + dZ^) , (3.27) 
a comparison with Eq. p.ip immediately reveals that 



B = ^f{R)h{R) , (3.28) 

w = 0, (3.29) 

V = ilog[/(i?)], (3.30) 

^i = Uog[h{R)]. (3.31) 



Inserting Eqs. ([5:^ - (P3T|) in Eqs. ^^-^?^, we find 

^ 2mR^J{R)h{Rf^{l + ^)'^ 

[VQU'iR) h{R) - f{R) h'm - f{R) h'{R)] , 

(3.32) 

where 

Q = l + 7P^ (3.33) 
h 

L = X X P, (3.34) 
in agreement with Eq. (5.7) in Ref. [23j . 



Let us now investigate how the Hamiltonian (|3.17l) is 
affected by a change of the radial coordinate R. Denoting 
the new radial coordinate by r = |a;| and defining 

J-'.§. (3,35) 

the radial derivatives of the metric potentials can be re- 
expresscd as 

fR = frJ, (3.36) 

where again / = B, uj, v, fi. The spin S, the derivatives 
of the metric potentials with respect to cos 6, and the 
quantities 

N = n = - , (3.37) 

r 

I = ez X N = X (3.38) 

V = v = nx^ (3.39) 

are not affected by the coordinate change. The same 
applies to the quantities P ■ V R and P ■ ^R appearing 
in Eqs. ((XT5)) . and ((X^ . In fact, in spherical 

coordinates, we have P-VR = —Pe sin 6 and P-^R = P^, 
hence 

PVR = pvr, (3.40) 
P-^R = p-^r, (3.41) 

where p ~ p/m and p is the conjugate momentum in the 
new coordinate system, i.e., pi = dX^ /dx^ Pj. On the 
contrary, since P ■ N = Pn, we have 

P-N^{p-n)J. (3.42) 

It is therefore straightforward to compute -ffg in a coor- 
dinate system related to quasi-isotropic coordinates by a 
rcscaling of the radius. Wc have 

Hs = Hso + Hss , (3.43) 

where 



Hi 



so 



WToWe ^ S2 (VQ + i) VQRH'' ^ 

+R{S ■ i-ir ip-vr) ( + 1 ) - ^J■cose{p • ") - ^/Q (i^r {p-vr) + {ficosO - Vcos0){p ■ 'n)i'^) 



Scosee^+'^(p -10(70 + 1) {S-n)e 



B' 



"{p-ir) (2^0 + 1) \.J vr R{S ■ v) - v^oso {S ■ n) e\B - J Br e^+'' {p ■ ir) (^Q + l) R{S-v) 



(3.44) 
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{p-vr){p-£,r){S ■£,)B + e^'^^'+''^ (p • ^ rf{S ■ v) 



2B {^/Q + \) VQRe \ 
+6^^ (i + Vq) VQ {S ■ v) + J {p ■ n) R[{p ■ vr){S ■ n) - J {p ■ n) R{S ■ v)]^^ B^ 

^2(p+,.) . ^ ^s-n) + e''+'^ J{p-n){p-^r)R{S-^)B 



e ^ Wcos e 



2B{^+l)^R-^\ 

{S ■ n) {p ■ vrf - J {p ■ n) R{S ■ v) {p ■ vr) - e^f" + ^ R^S ■ n) bA , (3.45) 



g = 1 + 7^' = 1 + e-2M (p . „)2 J2 



52 ^2 ^2 ' 



(3.46) 



and where R must of course be expressed in terms of the 
new radial coordinate r. 



IV. HAMILTONIAN FOR A SPINNING 
TEST-PARTICLE IN KERR SPACETIME IN 
BOYER-LINDQUIST COORDINATES 

In this section, we will specialize the Haniiltonian de- 
rived in the previous section to the case of Kerr spacetime 
in Boyer-Lindquist coordinates. 

We start from the metric potentials appearing in 
Eq. (|3.ip . which in the case of a Kerr spacetime take 
the form El 



B = 



o2f 



.2/i 



R ' 
2aMr 

A 
AS 

E 

R^' 



with 



A = + - 2Mr, 
2 



A 



2 I 2 
r + a , 

vo^ - a^Asin^ 



(4.1) 
(4.2) 
(4.3) 
(4.4) 

(4.5) 
(4.6) 
(4.7) 
(4.8) 



where the parameter a, which has the dimensions of a 
length, is related to the spin vector S'Kerr of the Kerr 
black hole by 



IS'kcitI 

M 



(4.9) 



The Boyer-Lindquist coordinate r is related to the quasi- 
isotropic coordinate R by 



R + M 



Rh 
R ' 



(4.10) 



where Rn ~ %/ AP — is the horizon's radius in quasi- 
isotropic coordinates. Note that the inverse of this trans- 
formation is given, outside the horizon, by 



R=^{r- M + y/A) 



(4.11) 



We then obtain that the derivatives of the metric poten- 
tials take the form 



B^ 



r - M 



A 



i?yA ' 

2aM[T.w^ - 2r2 (S + tn^)] 



A2 



r- M r 



2r vd"^ ~ {r ~ M) srn^ 6 



A 



fir 
BcosB 

^cose 



Mcos e — 
and we also have 



r 

0, 



4a^ Mr A cos 6 
~(AE-f 2Afrw2)2 

2a'^Mrvj^ cos 9 
(AS-f 2Mrtn2)S ' 
cos 6 



dR R 



dr 



(4.12a) 
(4.12b) 

(4.12c) 

(4.12d) 
(4.12e) 
(4.12f) 

(4.12g) 

(4.12h) 

(4.13) 



Inserting Eqs. (|4.12ap - (|4.12hp and Eq. (|4TT3l) into 
Eqs. (|3.43p - p.46p . we find that R cancels out both in 
Q, that is 



_ ^ ^ A{p-n)^ ^ (j3-$r)^S ^ (p-vr) 



A sin^ ( 



E sin^ e 



, (4.14) 
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and in the Hamiltonian Hs- In conclusion, the Hamil- 
tonian of a spinning test-particle in Kerr spacetinie in 
Boyer-Lindquist coordinates is 

H^Ht,s + Hs, (4.15) 

with 



Hns = P'pi + aJm^ + 7^' p,pj , (4.16) 



where a, f3^ and Y'' ^^e given in Eqs. (|3.6p - (|3.8p and 
need to be computed using the Kerr metric coefficients 
and with 

Hs^Hso+Hss, (4.17) 

where 



Hso — 



(^/Q+l) VQe 



1 - A'cose {p-n)£. 



'Q {^r ip-vr) + (/^cose - Vcose) {p ■ n) ^^) 
-JBref'+'^ip-^r) {S-v)\ , 



B' 



B 



(4.18) 



and 



^^SS — W (5 • 5Kcrr) 



- e~^+'' {p-vr){p-^r){S-^)B + e^^^'+'''> {p ■ $ rf [S ■ v) 



2B (VQ + 1) VQe , 

fl + Vo) \fQ{S-v)eB^ + J{P-n) \{p-vr){S-n)-J{p-n){S-v) 



e ^ LOcos e 



(p ■ I r)2 {S-n) + 6^+'' J (p • n) (p • I r) (S • 1)5 



25 (VQ + 1) ^/0 [ 

+ [iS ■ n) (p ■ vrf - J {p ■ n) (S ■ v) {p ■ vr) - e^i"- (l + v^) ^ {S ■ n) ^^'^ B^j ^ 

I 



(4.19) 



where we define 



B = 


BR = 




(4.20) 






r - il/ - 


(4.21) 


Br = 


Bj. R = 










(4.22) 


J = 


Ji? = ^ 




(4.23) 


KCIT = 


•S'kcit 




(4.24) 



nian for a spinning test-particle in a Schwarzschild space- 
time in Schwarzschild spherical coordinates by setting 
a = 0, and c omp are the result to the expression com- 
puted in Ref. Q [see Eq. (5.12) therein]. We find 



where 



i?3V0(i + VQ) 

M 



and we recall that — siv? 0. Wc stress that because 
this Hamiltonian is expressed in terms of quantities which 
are scalar under spatial rotations, we can express it in 
a cartesian coordinate system in which the spin of the 
Kerr black hole is not directed along the z-axis. For 
that purpose, it is sufficient to replace r with {x^ + 2/^ + 
z^)^/^, COS0 with ^Kerr ' 6^ with S'koit in Eq. (|3.38p . 
and express the vectors appearing in Eqs. (|4.16p - (j4.19p 
in terms of their cartesian components. 

As a consistency check, wc can compute the Hamilto- 



1 - 
S* 

R 



2R 



2 1- 



M 



{L-S* 



m 



M 
2R 



(4.25) 

(4.26) 
(4.27) 
(4.28) 



and 



r J sm 9 
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in agreement with Ref. [231. Also, it is worth noting 
that the Hamiltonian (|4.25p is the same as the quasi- 
isotropic Schwarzschild Hamiltonian ()3.32|) , expressed in 
terms of the Schwarzschild coordinate r. This is because 
the scalar product L ■ S \s unaffected by a change of the 
radial coordinate. 



descriptions, and write the improved EOB Hamil- 
tonian (see Sec. IV Ep . 



The ADM Hamiltonian canonically transformed 
to EOB coordinates 



EFFECTIVE-ONE-BODY HAMILTONIAN 
FOR TWO SPINNING BLACK HOLES 



We denote the ADM canonical variables in the binary's 
center-of-mass frame with r' and p' . It is convenient to 
introduce the following spin combinations: 



The EOB approach was originally introduced in 
Refs. [iO, [Hi to provide us with an improved (re- 
summed) Hamiltonian that could be used to evolve a 
binary system not only during the long inspiral, but also 
during the plunge, and that could supply a natural mo- 
ment at which to switch from the two body description 
to the one-body description, in which the system is rep- 
resented by a superposition of quasi-normal modes of the 
remnant black hole. 

A crucial ingredient of the EOB approach is the real 
PN-expanded Arnowitt-Deser-Misner (ADM) Hamilto- 
nian (or real Hamiltonian) describing two black holes of 
masses mi,m2 and spins S^i, ^2. The real Hamiltonian 
is then canonically transformed and subsequently mapped 
to an effective Hamiltonian H^s describing a test-particle 
of mass /i ~ mi 7712/ {mi + TO2) and suitable spin S*, 
moving in a deformed Kerr metric of mass M = mi +m2 
and suitable spin S'Kcn- The parameter regulating the 
deformation is the symmetric mass ratio of the binary, 
77 = fi/M, which ensures that the deformation disappears 
in the case of extreme mass-ratio binaries. The resulting 
improved EOB Hamiltonian then takes the form 



H- 



improved 



real 




(5.1) 



The computation of the improved EOB Hamiltonian 
consists of several stages. For this reason, we briefly re- 
view here the main steps and the underpinning logic that 
we will follow in the rest of this section: 

(i) We apply a canonical transformation to the PN- 
expanded ADM Hamiltonian using a generating 
function which is compatible with the one used in 
previous EOB work, obtaining the PN-expanded 
Hamiltonian in EOB canonical coordinates (see 
Sec. rvX)) : 

(ii) We compute the effective Hamiltonian correspond- 
ing to the canonically transformed PN-expanded 
ADM Hamiltonian (see Sec. IVB|) : 

(iii) We deform the Hamiltonian of a spinning test- 
particle in Kerr derived in Sec. IIVI by deforming 
the Kerr metric (see Sec. IV Cp , and expand this 
deformed Hamiltonian in PN orders (see Sec. IV Dl : 

(iv) Comparing (iii) and (iii) , we work out the mapping 
between the spin variables in the real and effective 



(To = (T 



S1+S2., 

ma 
cr*. 



^1 



7712 



(5.2) 
(5.3) 
(5.4) 



Moreover, in order to consistently keep track of the PN 
orders, we will restore the speed of light c and rescale the 
spins variables as cr* — > cr* c and cr ^ cr c.^ The canoni- 
cal ADM Hamiltonian is known through 3PN order [26l - 
l30j and partially at higher PN orders [3l|, HI]. In par- 
ticular, the spin-orbit and spin-spin coupling terms agree 
with those computed via effective-field-theory techniques 
at 1.5PN, 2PN and 3PN order In this paper, we 

use the spin-independent part of the ADM Hamiltonian 
through 3PN order, but we only use its spin-dependent 
part through 2.5 PN order, i.e., we consider the leading- 
order (1.5 PN) and the next-to-leading order (2.5PN) 
spin-orbit couplings, but only the leading order (2PN) 
spin-spin coupling. The expressions for these couplings 
are [Illl^ 

1 T' 

Hi§^ir',p',a*,a) ^ ^ ±^ . (g^^^ ^ + g^P^ a*) , 

(5.5) 

i?4°^(r',p',^*,<T) ^ [3{n' ■ aof - al] , 



(5.6) 



with L' = r' X p', n' = r'/r', and 



ADM 

9(7 



^ADM 



2 + ^ 



2 



^rjp'^ + ^ri{n' -p'f 



(5.7a) 



3 1 
2 ? 



^ \ -/2 I 3 , , ,,2 



-(5 + 27?)- 



(5.7b) 



^ This is appropriate for black holes or a rapidly rotating com- 
pact stars. In the black-hole case, S = X-'^^^/c, with x rang- 
ing from to 1. In the rapidly spinning star case one has 
S = Mvrotr ~ McTs ~ M^/c (where we have assumed that the 
rotational velocity Vrot is comparable to c and that the stellar 
radius r is of the order of the Schwarzschild radius Ts ~ M/c^). 
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where we have introduced the rescaled conjugate momen- 
tum p' = p' / ^i. 

We now perform a canonical transformation from the 
ADM canonical variables r' and p' to the EOB canonical 
variables r and p. Let us first consider the purely orbital 
generating function 

G(r',p)=r'-p + GNs(r-',p), (5.8) 

GNs(r',p) =GNSlPN(r',p) 

+ Gns2Pn(^',p) +Gns3Pn(^',p) , (5.9) 



where the IPN-accurate generating function Gns ipn was 
derived in Ref. [l^l , 



G 



NSIPN 



(r',p) 



1 



-r p 



1 



M 



1 



1 



(5.10) 

while the 2PN and 3PN accurate generatiiie functions, 
Gns2PN and GnsspN: were derived in Refs. [To| and 
respectively. From the definition of generating function, 
it follows that the transformation of the phase-space vari- 
ables is implicitly given by 



Pi = Pi 



dG^s{x',p) 
dpi 

dx"- 



(5.11) 
(5.12) 



while the Hamiltonian transforms as H{r,p) = 
H^^^{r' ,p'). At linear order, which is enough for 
our purposes, Eqs. (|5.1ip and (|5.12p can be written as 
y = y' — {Gns, 2/'}, where {...} are the Poisson brackets 
and where y stands for either x or p. The transforma- 
tion of the Hamiltonian, again at linear order, is then 
H{y) = H^^^{y) + {Gns, i/^°^}(y) 0- Similarly, if 
one considers a generating function which depends not 
only on the orbital variables, but also on the spins, 

G(r',p, (T* ,(t) = r' ■ p 

+ GNs(r',p) + Gs(r',p,rT*,^) (5.13) 

the Hamiltonian will again transform as H{y) ~ 
^ADM(y) ^ {GNs,i?^°^}(y) + {Gs,i/AD^}(2/), where 
now the Poisson brackets in the term {Gs, H^^^} will 
involve also the spin variables [2l|. In particular, let us 
consider a spin-dependent generating function 

Gs{r',p,(T*,eT) =Gs2PN(T-',p,cr) 

+ Gs2.5PN(r',p,cr*,(T) (5.14) 
+ Gsss2.5Pn(t-',p,o-*,(t) . 

where the 2PN-accuratc spin-dependent generating func- 



tion Gs 2PN was implicitly^ used in Ref. [l9[ , 
Gs2PN(r',p,<T) = - {[a' (a ■ n')^](r' -p) 



+ (cr • n'){r' x p) ■ {cr x n')| 



(5.15) 



the 2.5PN-accurate generating function Gs 2.5PN linear in 
the spin variables was introduced in Ref. |2l| . 



Gs2.5PN(r',p,cr*,(T) 



(r' • p){r' X p) 



Kt?) <T + &(77) <T*] , (5.16) 



a{rf) and h{if} being arbitrary gauge functions; also, for 
reasons which will become clear in Sec. IV Dl we include 
the following 2.5PN-accurate generating function, cubic 
in the spins. 



GsSS2.5PN(T"',P,cr*,'7') = 



2Af3r'4c5 



{(T-r')[cT* -{cTXr')] 
(5.17) 



When applying the generating function (|5.13p to the 
ADM 2PN spin-spin Hamiltonian (j5.6|) . we obtain 



i^ss2PN(r,p,<T*,^) = H^i>^^{r,p,cT\<r) 

+ {Gs2PN,-ffNcwt}(^,P,0-) 



with 



Newt — 



r 2/z 



(5.18) 



(5.19) 



{Gs2PN,i?Ncwt}(r,p,cr) = — - — [{n ■ erf - a- 



+ [{p-2{p-n)n)-(T]p-(TY (5.20) 

Similarly, if we apply the same generating function to 
the ADM spin-orbit Hamiltonian (|5.5p . the 1.5PN order 
term remains unaltered [2l|, while the 2.5PN order term 
transforms as [2l] | 

^^S0 2.5PN {r,p,(T*,cr) = H^§'^^p^{r,p,cr*,cr) 
+ {G2.5PN,-ffNcwt}(r,p,cr*,cr) , 

+ {Gns ipn, iJg'oTsPNK^, P, (5-21) 

where 



G 



2.5PN 



G; 



SS2.5PN 



G 



SSS2.5PN 



(5.22) 



See discussion in Sec II D of Ref. The need for this gener- 

ating function will become apparent with Eq. I I5.55I I in Sec. IV E)] 
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{G2.5PN, H^owt}ir,p, (T*, (t) 



1 



L-[b{7i)rT*+a{7^)rT] 



hp - 3(p • n) 



[cr* • (cr X n)] [cr ■ [p - 2{p ■ n) n)] 
(i • cr*)cr2 {L ■ (t){(t* ■ cr) 



(5.23) 



and 



{GnS IPN, ^So'i^SPnK'": Pt^*i^) 

3L /3 



(T* + 2^^ |-^(2 + j;)+7?[p2 + 2(p-n)2]| 

(5.24) 



Therefore, the complete real Hamiltonian in the EOB 
canonical coordinates is 

H{r,p,(T* ,(t) = i?nospin(r',P, cr*,cr) 
+H^S''{r,p,cT\cr) 
+Hii^^{r,p,cT*,cT) 

+ {G2.5PN, i?Nowt}(r, p, (T*,(t) 

+{Gnsipn, Hs§^r,p^}{r,p, cr*, cr) 

+ {Gs2PN,-ffNowt}(T',P,'T) , 

(5.25) 

where iJnospin is the 3PN ADM Hamiltonian for non- 
spinning black holes, canonically transformed to EOB 
coordinates, which can be obtained from Ref. [TT| . 



B. Spin couplings in the effective Hamiltonian 

Following Refs. [111, [l^, we map the effective and 
real two-body Hamiltonians as 



2m 1 m2 



(5.26) 



where -ffjcai is the real two-body Hamiltonian containing 
also the rest-mass contribution M (? . We denote the non- 
relativistic part of the real Hamiltonian by H^^, i.e., 
H^ ^ EE ff i-eai-M c^. Identifyi ng H^ ^ with H as given in 
Eq. (|5.25p , and expanding Eq. (|5.26p in powers of 1/c, we 
find that the 1.5PN and 2.5PN order spin-orbit couplings 
of the effective Hamiltonian are 



H^^{r,p,a*,rr) 



1 L 



[cr* • (cr X n)] [cr ■ {p — 2{p ■ n) n)] 
{L ■ a*)a^ - {L ■ ct){(7* • cr) 

(5.27) 



where (211 



9a 



2 H 



-1] + a{ri) 

{p-nf 



-77 + 3 a(77) 



(5.28a) 



1 



1 



8 



77 + 6(77) 



15 



77 + 3 6(77) 



[p-nf 



M 
r 



■ 7] + b{r]) 



2 4 

and the 2PN order spin-spin coupling is 



(5.28b) 



1 7; 



2r 



[{p-2{p-n)n)-(T]p-a-Y (5.29) 



C. The Hamiltonian of a spinning test-particle in a 
deformed Kerr spacetime 

We now deform the Hamiltonian of a spinning test- 
particle in a Kerr spacetime computed in Sec. IIVI [see 
Eqs. (gm, (|iT7l) . and (|iT^ ] by deforming the 

Kerr metric. The deformation that we introduce is reg- 
ulated by the parameter rj = /x/M, and therefore dis- 
appears in the test-particle limit. Also, the deformed 
Hamiltonian will be such as to reproduce, when expanded 
in PN orders, the spin couplings of the effective Hamil- 
tonian given in Sec. IV Bl 

When the spin of the Kerr black hole is zero, that is a = 
0, we require the metric to coincide with the deformed- 
Schwarzschild metric used in the EOB formalism for non- 
spinning black- hole binaries [13, [13 . That deformation 
simply amounts to changing the components gu and ff^r 
of the metric. In the spinning case, following Ref. [21|, 
we seek an extension of this deformation by changing the 
potential A appearing in the Kerr potentials (|4.1l) - (|4.4p . 

It is worth noting, however, that we are not allowed 
to deform the Kerr metric in an arbitrary way. We re- 
call indeed that the Hamiltonian that we have derived 
in Sec. IIVI is only valid for a stationary axisymmetric 
metric, and in coordinates which are related to quasi- 
isotropic coordinates by a redefinition of the radius. In 
other words, it must be possible for our deformed metric 
to be put in the form (|3.ip by a coordinate change of the 
type R = R{f)- For this reason we cannot deform the 
metric exactly in the same way as in Ref. |2l| . Here we 
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propose to deform the metric potentials in the following 
manner 



B 



UJ 



o2m _ 



Ji 

At ' 
At 

^2 ' 



and 



J-i = — = 
dr 



R 



(5.30) 
(5.31) 
(5.32) 
(5.33) 

(5.34) 



where the relation between r and i? can be found by 
integrating Eq. (|5.34p : 



R = exp 



dr 



The deformed metric therefore takes the form 



5" 


A 












1 










9** 





At E sin^ e 
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t<t> 



"AtE 



(5.35) 

(5.36a) 
(5.36b) 
(5.36c) 
(5.36d) 

(5.36e) 



which does not depend on R. Therefore, as we will show 
explicitly later in this section, we do not need to com- 
pute the integral (|5.35p to write the Hamiltonian. The 
quantities At, A^, At and uifd in Eqs. (|5.36ap - ()5.36ep are 
given by 



At = ^A{u) + ^u' 

Ar - AtD-\u), 

At = tu'* - At sin^ , 



2a M r + ?/ 



(5.37) 

(5.38) 
(5.39) 

(5.40) 



where u = M/r, oj^^ and uj^^ are adjustable parameters 
which regulate the st reng th of the frame-dragging, and 
through 3PN order d, [Tj 



A(?i) = l-2u + 2r]u^ + 7] i^— - —T 
D-^(u) = l-|-6?/u2-|-2(26-37?)?7M^. 



94 41 2> 4 



(5.41) 
(5.42) 



We find that our deformed metric is the same as the de- 
formed metric of Ref. [2l|, except for g"^"^ and 5*"^.^ As 
we prove below, the differences between our deformation 
and the deformation of Ref. [2l| appear in the Hamilto- 
nian at PN orders higher than 3PN. 

To obtain the total Hamiltonian ()4.15p . that is H = 
^NS + I we first compute the Hamiltonian -Sns for a- 
non-spinning particle in the deformed-Kerr metric. Using 
Eq. (|4.16p and Ref. [ll|, we have 



/3' p, + a^m^ + 7^' p^ Pj + Qi{p) , (5.43) 



where Qiij)) is a term which is quartic in the space mo- 
menta Pi and which was introduced in Ref. [ll| , and 



a 



V- 



9'' 9'' 



(5.44) 

(5.45) 
(5.46) 



In Eqs. (|5.44p - (|5.46p the metric components have to be 
replaced with those of the deformed-Kerr metric (|5.36ap - 
(|5.36ep . When expanded in PN orders, Eq. (|5.43p co- 
incides, through 3PN order, with the Hamiltonian of a 
non-spinning test particle in the deformed-Kerr metric 



given by Ref. [21 



Second, to calculate iJs given by Eqs. (jXi^l) . ([5^ 
and p.45p . we need to compute the derivatives of the 
metric potentials. We obtain 



Br 



fA'rAt' -2At 



2^/K;rK'tR 

-K[ wfd -I- At 



A? 



Mr 



Br, 



r 
S 
r 

0, 



'At'-4rAt 



2At At 



1 



Wcos e — 



IJ-cosB 



/Ar 

2a} cos B At uj^a 



A? 



a w cos d{vj^ — At) 



AtE 



a} cos( 



(5.47a) 
(5.47b) 

(5.47c) 

(5.47d) 
(5.47e) 
(5.47f) 

(5.47g) 
(5.47h) 



where the prime denotes derivatives with respect to r. 
As already stressed, although the metric potentials _B, 



Ref. [211 chooses 



1*0 



(-a^ sin^e -I- At)/(AiSsin2 6») and 



a (At — ro )/(At S), which are different from our expres- 



sions I l5.36dl l and I l5.36el l even for tj! 



,fd 



: 0. 
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w, V and /i depend on R, the factors i? cancel out in 
the deformed-Kerr metric. Therefore, those factors must 
cancel out also in H^. This happens because the ref- 
erence tetrad field e-A which, together with the metric, 
completely determines the Hamiltonian [see Eq. p.9p ]. 
can be defined independently of R. Indeed, this turns 
out to be the case, and if we introduce the rescaled po- 
tentials 



B 

Br 

J 




(5.48) 

(5.49) 

(5.50) 
(5.51) 



the case of the Kerr metric, and they are given by 
1 2 



NS 

SO 1.5PN 



L ■ Skc 



NS 

S0 2.5PN 



rrNS 
^SS2PN 



(5.53) 

(5.54) 



^3 j,3 

= 0, 
1 in 

+ [{p-2{p- n)n) ■ Skcit] P ■ S^Kcrrj • 

(5.55) 



and define 



Q = 1 



A.(p-n)2 (p.|r)2E 



(P 



V r)^ 



At sin^ S sin^ 6 



(5.52) 



the Hamiltonian Hg for the deformed-Kerr metric takes 
exactly the same form as in the Kerr case [see Eqs. (|4.17p . 
(|4T8l) and (|4?T9| ]. where we recall that = sin^ and 
where now uj and its derivatives, and its derivatives, and 
the derivatives of /i are given by Eqs. (|5.3ip . (|5.32p . and 
Eqs. ([5:47a|) - (|5.47h|) . Also, as we have already stressed, 
in order to express the Hamiltonian in a cartesian 
coordinate system in which the spin of the deformed-Kerr 
black hole is not directed along the z-axis, it is sufficient 
to replace r with (a;^ + + z'^Y^^, cos6 with SKcn- ■ n, 
Bz with (^Kerr in Eq. (|3.38p . and to express the vectors 
appearing in the Hamiltonian in terms of their cartesian 
components. 



D. PN expansion of the deformed Hamiltonian 

We now expand the deformed Hamiltonian H = Hjqg + 
Hs derived in the previous section into PN orders. We 
will denote the spin of the deformed-Kerr metric with 
S^Kcrr, while for the test particle's spin we introduce the 
rescaled spin vector S* = S M/m, S being the physical, 
unrescaled spin. Also, we rescale the spins as SKcn ^ 
5^Kcrr c and S* S* c, so as to keep track of the PN 
orders correctly. Moreover, we set S'xerr = XKerrM^, 
XKcrr being the dimensionless spin of the deformed-Kerr 
black hole, with norm |xKcrr| ranging from to 1. 

As already mentioned, the part of the Hamiltonian 
which does not depend on the test particle's spin, HnSi 
agrees through 3PN order with the corresponding H^s 
computed in Ref. [2l|. Moreover, although the metric 
(|5.36ap - (|5.36ep only coincides with the Kerr metric for 
77 = 0, the dependence on rj appears neither in the 2PN 
order coupling of the deformed-Kerr black hole's spin 
with itself, nor in its 1.5PN and 2.5PN order spin-orbit 
couplings. Those couplings are therefore the same as in 



Expanding then in PN orders the part of the Hamiltonian 
that depends on the test particle's spin, that is Hg, we 
find 



^SO 1.5PN 



FfS 

-"S02.5PN 



2r3 c3 

1 



L-S* , 

M fl 



(5.56) 



L-S* 



[S* ■ (^Kcrr X w)] [^Kcrr ' {p - 2{p ■ u) u)] 
{L ■ S*)S^^^.j. — {L ■ ^Kcrr) {S* ■ Skcit) 



Hi 



SS2PN 



Mr3 



2M3 r4 c5 



{Srii Uj Sij ) S'kcit 



(5.57) 
(5.58) 



We recall that the Hamiltonian for a spinning test parti- 
cle in curved spacetime from which we started the deriva- 
tion of our novel EOB model [see Eq. p.9p ] is only valid 
at linear order in the particle's spin. Therefore, the 
same restriction applies to the Hamiltonian derived in 
Sec. IV CI In particular, that Hamiltonian docs not in- 
clude the couplings of the particle's spin with itself. We 
introduce those couplings by hand, at least at the leading 
order (2PN), by adding a quadrupole deformation [l^l 
/i^'', quadratic in the particle's spin, to the deformed- 
Kerr metric in Sec. IVH] [see Eqs. ([05a|) - (iOSc)) ]. The 
expression for /i'*" and the details of the above proce- 
dure — together with a way in which it can in principle 
be extended to reproduce also the next-to-leading order 
coupling of the particle's spin with itself — are given 
in Appendix [X] For the purpose of the present discus- 
sion, however, it is sufficient to mention that the addition 
of this quadrupole deformation to the metric (|5.36ap - 
(|5.36cp augments Eq. (|5.55p by the term 



-(3nj Uj ~dij)Sl Si 



(5.59) 



2M r3 c4 

Therefore, the total leading order spin-spin Hamiltonian 



12 



E. The efFective-one-body Hamiltonian 



HsS 2PN — H^s 2PN + H^S2PN 



2Mr3 c4 
m 

2Mr3 c4 
1 m 



(3n, rij -(5,j)5* 5^ 



[(p - 2{p ■ n)n) ■ ^Kcrr] P ■ SKcrr} , (5.60) 



With 5^ = 



As we will show in Sec. IV El a proper choice of the 
vectors S^Kcn- and S* in terms of the vectors cr and cr*, 
defined in Eqs. (|5.2p and (|5.3p . allows us to reproduce 
the PN-expanded effective Hamiltonian [see Eqs. (|5.27p - 
(|5.29p ] using the PN-cxpandcd dcformcd-Kcrr Hamilto- 
nian that we have just derived. 

Finally, it is worth noting that the presence of terms 
quadratic in the deformcd-Kcrr black hole's spin in 
Eq. (|5.57p explains why we introduced the 2.5PN- 
accurate canonical transformation (|5.17p . Indeed, the 
latter produces exactly the same terms in the PN- 
expanded effective Hamiltonian (|5.27p at 2.5PN order. 
Quite interestingly, the terms quadratic in S'kcit appear- 
ing in Eq. (|5.57p could also be eliminated with a suitable 
choice of the reference tetrad ca- In fact, as stressed in 
Sec. mil and in Ref. a choice of the reference tetrad 
field corresponds to choosing a particular gauge for the 
particle's spin. In agreement with this interpretation, 
we find that the terms of Eq. (|5.57p which are quadratic 
in ^Korr disappear if the initial tetrad (|3.16ap - (j3.16dp is 
changed to a different tetrad related to the original 
one by the following purely-spatial rotation: 



where the rotation matrix TZij is given by 



n = nY 



a^XZ 



2i?4 



Tlx 



■Y Z 



2i?4 



(5.61) 



(5.62) 



TZxii^] and 7?.y[(/)] being rotations of angles ip and 4> 
around the axis X and Y, respectively. 

As a consistency check, we have verified that this new 
tetrad is the same as that used in Ref. [2^ when com- 
puting the Hamiltonian in ADM coordinates, where those 
terms quadratic in S^Kcn- do not appear. We have checked 
this by transforming the new tetrad (|5.6ip from quasi- 
isotropic to ADM coordinates [which are related by the 
coordinate transformation (49) in Ref. [3l|], and compar- 
ing it to the tetrad given in Eqs. (6.9a)-(6.9b) of Ref. |23| . 
and find that the two tetrads agree through order 1 /c^ . 



In this section we first find the mapping between 
the masses /i, M and the spins <t and cr* of the ef- 
fective Hamiltonian derived in Sec IV Bl and those of 
the deformed-Kerr Hamiltonian derived in Sees. IV CI 
and IV Dl that is m, M, S'Korr and S* . Then, we derive 
the improved (resummed) EOB Hamiltonian. 

As shown in Ref. Q , matching the non-spinning parts 
H^s of these Hamiltonians forces us to identify the total 
mass AI of the two black holes in the PN description with 
the deformed-Kerr mass M of the test-particle descrip- 
tion, thus justifying our choice of using the same symbol 
for these two a priori distinct quantities. Similarly, we 
find that m = /i [9| . Assuming this mapping between the 
masses and imposing that the PN-expanded deformed- 
Kerr Hamiltonian given by Eqs. (|5.53p - (|5.60p coincides 
with the effective Hamiltonian given by Eqs. (|5.27p - 
(|5.29p . we obtain the following mapping between the 
spins 



S* 



(5.63) 
(5.64) 



where we have set for simplicity a(?/) — and 6(77) = 
and where 



+6(p- 



^(4.-7.*) 
r 



■5cr*) 



p2 (3^. 



4cr* 



(5.65) 



Here, Ao-- is an arbitrary function going to zero at least 
linearly in rj when 77 — > 0, so as to get the correct 
test-particle limit. In fact, if Ag-* satisfies this condi- 
tion and if we assume, as appropriate for black holes, 
5^1,2 = Xi,2'iTT'i2 (with |xi,2| < 1 and constant)^, when 
TO2 ^ wc have Skcti- = Si + 0{m2). Similarly, for 
m2 the physical unrescaled spin of the effective par- 
ticle is S = S* m/M ^82 + 0{m2 f. The equations of 
motion of our initial Hamiltonian p.9p coincide with the 
Papapetrou equations [2^ , which describe the motion of 
a spinning test-particle in a curved spacetime [s^, [3^ . 
Assuming the canonical commutation relations between 
x% Pj, S'l and S'2, we obtain that the Hamilton equa- 
tions for the effective deformed-Kerr Hamiltonian are 
y ~ yp + 0{m2). Here, the dot denotes a time deriva- 
tive, y is a generic phase-space variable (a;% pj, Si or 



® As noted by Ref. [2l| . a spin mapping such as ours also gives 
the correct test particle limit if |Si,2|/'"i,2 = const., but this 
scalin g of the spins with the masses is not appropriate for black 
holes 
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S'2), and y = yp are the Papapetrou equations expressed 
in Hamiltonian form. Therefore, our mapping repro- 
duces the correct test-particle hmit, and the remainders 
'S'Korr— S'l = 0(7712) and S—S2 = 0{m2f' produce extra- 
accelerations of order 0(7772) or higher. This is compa- 
rable to the self- force acceleration [40| , which appears at 
the next order in the mass ratio beyond the test-particle 
limit. 

Although different choices for the function Ag-* arc in 
principle possible, we choose here 



A.. = ^ 



r 



- p2 (3cr 



4cr*) 



-6(p-n)2(6cr + 5cr*) 



(5.66) 



which gives, when inserted into Eq. (|5.65p . A(j = 0. Be- 
cause this form for Ao-* is clearly not covariant under 
generic coordinate transformations, we choose instead 
the following form for the mapping of the spins, which is 
covariant at least as far as the square of the momentum 
is concerned: 



0, 



A„. 



, , 
{In* 

r 



-6-(p 



^ 4(t) + {Q-l) (3cr - 
nf (6cr + 5cr*)j , 



where we have replaced with j^-'pi'Pj = Q 
Eq. (ESg)] and {p ■ nf = 



(5.67) 

4(T*) 

(5.68) 
1 [where 

Q is given in Eq. (|5.52p ] and {p ■ nY = with 
Ar(p • rif' /Tj = g^^pf.. This form agrees with the pre- 
vious mapping through order 1 /c^ , but differs from it at 
higher orders. Although neither this form is completely 
covariant, not even under a rescaling of the radial coor- 
dinate (as it still features a dependence on the radius 7-), 
it proved slightly better as far as the dynamics of the 
EOB model, analyzed in the next section, is concerned. 
In particular, the factor g^^, which becomes zero at the 
horizon, quenches the increase of Pr at small radii, thus 
giving a more stable behavior during the plunge subse- 
quent to the inspiral. (A similar effect was observed in 
Rcf. [2^, where the radial momentum was expressed in 
tortoise coordinates to prevent it from diverging close to 
the horizon.) 

Having determined the mass and spin mappings, we 
can write down the improved (rcsummed) Hamiltonian 
(or EOB Hamiltonian) for spinning black holes. To this 
purpose, it is sufficient to invert the mapping between 
the real and effective Hamiltonians [Eq. (|5.26|) ]. In units 
in which c = 1, we obtain 



H. 



improved 



real 




(5.69) 



with 



Hs + 13" Pi + a ^Jn^+Y^PiPj + Q^ip) 
{S'^ ^3n'n^)S* S* . 



2Mr^ 



Here, the -n/{2M r^)(6^^ - 3n^n^)S*S* term is the 
quadrupole deformation introduced in the previous sec- 
tion to account for the leading order coupling of the par- 
ticle's spin with itself (see also Appendix IX)) : /?% a and 
7*-' are computed using the deformed-Kerr metric, that 
is inserting Eqs. (|5.36a|) - (j5.36cp into Eqs. 
Hs is obtained by inserting Eqs. (|5.3ip . (|5.32l) . and 
Eqs. (|5:T7a|) -(|g3 ^ into Eqs. (|iT7)) . and gH]). 

Lastly, the spin S'Kerr enters this Hamiltonian through 
the parameter a = \Skct:i-\/M appearing in the deformed- 
Kerr metric. 

Before completing this section, we want to discuss the 
deformation of the Kerr potentials At and A^ given in 
Eqs. (|5.37|) and (|5.38p . which play an important role in 
the EOB Hamiltonian (j5.69p . It is convenient to re- write 
the function At as 



A„(u) = A{u) + —, 



(5.71) 
(5.72) 



In previous EOB investigations the Fade summation was 
applied to the function to enforce the presence of a 
zero, corresponding to the EOB horizon, both in the non- 
spinning 11 1 and spinning case [l^ [21[ . Reference [22j 
pointed out that when including the 4PN and 5PN terms 
in the function A{u)^ the Padc summation generates 
poles if spins are present. Also, the Padc summation docs 
not always ensure the existence of an innermost stable 
circular orbit (ISCO) for spins aligned and antialigned 
with the orbital angular momentum and, even when it 
does, the position of the ISCO does not vary monotoni- 
cally with the magnitude of the spins. For these reasons, 
we propose here an alternative way of enforcing the exis- 
tence of the EOB horizons. Working through 3PN order, 
we write 



Au{u) = A„(u) [l + ?7Ao-f-77log(l + AiM-^A2; 



+A^u^ + A^u^)] , 



where 



Au{u) 



^EOB 
H,± 




M2 -qK - 1 {r)K - 1)2 ' 
(a/ ± -a^^ {l~K if) . 



(5.73) 

(5.74) 

(5.75) 
(5.76) 



(5.70) 



Here, are the EOB horizons, which differ from the 

Kerr horizons when the adjustable parameter K is differ- 
ent from zero, and where the log is introduced to quench 
the divergence of the powers of u at small radii. We 
could in principle replace the logarithm with any other 
analytical function with no zeros (e.g., an exponential). 
However, when studying the dynamics of the EOB model 
(see Sec. IVI[) the results are more sensible if we choose a 
function, such as the logarithm, which softens the diver- 
gence of the truncated PN series. 
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The coefficients Aq, Ai, A2, A3 and A4 can be de- and (j5.72p . with A{u) given by its PN expansion (|5.41l) . 
rived by inserting Eq. (|5.73p into Eq. (|5.7ip . expanding Doing so, we obtain 
through 3PN order, and equating the result to Eqs. ()5.7ip 



Ao 
Ai 

A2 
A3 
A4 



Ai (-4?/ A' + Ai + 4) 



-2{rjK -l){K + Ao), 
1 
2 
1 
3 
] 

12 



IVP 



- Af + 3(?7 A' - 1) Af + 3A2 Ai - 6(77 A' - 1) (-?? A' + A2 + 1) - 3— {rj K - ly A 



^ {e-^ (A^ - 2A2) (77 K - if + 3At -8{r^K- 1) A? - 12A2 A^ + 12 [2(77 - 1) A2 + A3] Ai 
12 (f - 1-^) iv K-lf + 6 [Al - 4A3 (77 if - 1)] } . 



(5.77) 
(5.78) 

(5.79) 
(5.80) 

(5.81) 



By construction, if we expand Eq. (|5.73p in PN orders, 
AT can only appear at 4PN and higher orders, because 
we must recover the PN expansion (|5.37p ~ (|5.4ip through 
3PN order. In this sense, K parameterizes our ignorance 
of the PN expansion at orders equal or higher than 4PN 
(i.e., K would not play any role if the PN series were 
known in its entirety). Similarly, we re- write the poten- 
tial A^ [Eq. ((08)) ] as 



A, = AtD-\u), 



(5.82) 



D-^{u) = l + log[l + 677 1*2 + 2(26 -377) 771*3]. 

(5.83) 

The coefficients in the above function D~^{u) are such 
that, when PN expanded, it gives the PN result (|5.42p . 
and the logarithmic dependence is once again chosen to 
quench the divergence of the truncated PN series. 

Finally, let us stress that if we included PN orders 
higher than 3PN in the functions A{u) and D{u), we 
would need to add higher order coefficients A^ with i > 4 
in Eq. (fSTTH]) . 



VI. EFFECTIVE-ONE-BODY DYNAMICS FOR 
CIRCULAR, EQUATORIAL ORBITS 

In this section we study the dynamics of the novel EOB 
model that we developed in Sec. IV El We will show that 

(i) Our EOB model has the correct tcst-particle limit, 
for both non-spinning and spinning black holes, for 
generic orbits and arbitrary spin orientations; 

(ii) There exist an IS CO when the spins are aligned or 
antialigned with the orbital angular momentum L; 

(iii) The radius, energy, total angular momentum, or- 
bital angular momentum and frequency at the 



ISCO exhibit a smooth dependence on the binary 
mass-ratio and spins. Also, this dependence looks 
reasonable based on what we expect from the test- 
particle limit and from numerical-relativity simula- 
tions; 

(iv) The frequency at the ISCO for an extreme mass- 
ratio non-spinning black-hole binary agrees with 



the exact result computed by Ref. |4l| : 



(v) During the plunge subsequent to the ISCO, the 
orbital frequency of black-hole binaries with spins 
aligned or antialigned with L grows and reaches 
a maximum, after which it decreases. The ra- 
dius at which the frequency peaks is very close to 
the radius of the equatorial, circular light ring (or 
photon orbit). This feature generalizes the non- 
spinning behavior and it has a clear physical 
interpretation in terms of frame-dragging. As in 
the non-spinning case , it provides a natural time 
at which to match the two-body description of the 
inspiral and plunge to the one-body description of 
the merger and ringdown. 

We stress that only (i) applies to generic orbits and spin 
orientations, while (ii), (iii), (iv) and (v) are true for 
black-hole binaries with spins aligned or antialigned with 
L. (It should be noted that circular or spherical orbits, 
and therefore the ISCO, arc not even present for generic 
orbits and spin orientations, because the system is not in- 
tegrable, not even in the test-particle limit [13]). While 
we will tackle the study of generic orbits and arbitrary 
spin orientations in a follow-up paper, we argue that the 
preliminary study presented here is already sufficient to 
illustrate the potential of the novel EOB model. We re- 
call that the only existing EOB model for spinning 
black-hole binaries, proposed in Refs. 0,1211, (0 repro- 
duces only approximately the test particle limit; (ii) when 
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including non-spinning terms at 4PN and 5PN order, it 
does not always present an IS CO for binaries with spins 
parallel to L, and when it does the spin dependence of 
quantities evaluated at the ISCO is unusual; (iii) gen- 
erally, the orbital frequency does not peak during the 
plunge, making the prediction of the matching time from 
the two-body to the one-body description quite problem- 
atic. 

Let us now go through the points of the list that we 
presented at the beginning of this section. In order to 
prove point (i) we first need to observe that the de- 
formed metric (|5.36ap - (|5.36ep [with the potentials 
and At given by Eqs. (|5.82|) and (|5.73|) ] reduces to the 
Kerr metric as 77 — >■ 0, and the deformation is linear in 
•q when 77 ~ 0. Therefore, the acceleration produced by 
this deformation on the test-particle is comparable to the 
self-force acceleration, which appears at the next order 
in the mass ratio beyond the test-particle limit. Sec- 
ond, as already proved in Sec. IV E| the mapping (|5.63|) - 



(|5.64[) of the spins reduces to S'xerr = Si 0(7712) and 
S = S*m/M =82 + 0(777,2)^ when 7712 0, where the 
remainders produce accelerations which are again com- 
parable to the self-force acceleration. 

To prove points (ii), (iii), (iv) and (iv), we need to 
write the effective EOB Hamiltonian (|5.70|) for equato- 
rial orbits and for spins parallel to the orbital angular 
momentum (chosen to be along the z-axis). We obtain 



M ^2 
2Mr3 * ' 

gfoL-S*+gfsS\ 



(6.1) 
(6.2) 



where 



J 



eff 

9so 



V^7^(B^ - 2B,^r) + {e^'+'' - (x/^+ l) + {Bu, - Br)y^ 



B2M(yQ + l) VQ 

- B^ e-2(A+.) A 



2B (VQ + i) VQ 



(6.3) 
(6.4) 



with 



^2 j,2 



^2 (tU^ 



l2 A, 



(6.5) 



The above equations can be evaluated explicitly by us- 
ing Eqs. dOl]), (IO2I) . (l5Tfi)) - (l5:47a . (lOSll - dtsOl) . 
(fSTTTI) - (fS^ "^. To calculate the radius and the or- 
bital angular momentum at the ISCO for the EOB 
model, we insert Eq. (|6.1|) into the real EOB Hamilto- 
nian (j5.69p . and solve numerically the following system 
of equations @ 



dH'::r^'^\r.Pr=Q.L^ 

dr 

do rrimprovcd / 



0, 



0,L.) 



dr 



0, 



(6.6) 
(6.7) 



with respect to r and L^. Moreover, the frequency for 
circular orbits is given by 



^rrimprovod/ 



0,iz) 



(6.8) 



A Mathcmatica notebook implementing the Hamiltonian II6.III — 
116.51 1 is available from the authors upon request. 



which follows immediately from the Hamilton equations 
because Lz = p^,- Finally the binding energy is Eund = 

improved 



H. 



real 



M. 



Henceforth, we set the adjustable frame-dragging pa- 
[sec Eq. (|5.40|) ] and write K in 



rameters u) 



fd 



,fd 



Eq. (j5.76p as a polynomial of second order in 7;, 
K{t]) = Ka + Kir] + K2 rj^ 



(6.9) 



Kq, Ki and K2 being constants. We find that if we 
impose 



K{l/4) = \ 







(6.10) 



the functional dependence on rj and x of several physical 
quantities evaluated at the ISCO is quite smooth and 
regular. Therefore, imposing these constraints we obtain 



KiT^) ^ Ko (1 - 477)2 + 4(1 - 2?/)?/ . 



(6.11) 



It is worth noting that the values of K and dK/drj at 77 = 
1/4 have a more direct meaning than the coefficients Ki 
and K2. In fact, current numerical-relativity simulations 
can evolve binary black holes with 77 « 0.25 (with only 
few runs having rj ~ 0.1). Thus, they can determine 
K{l/4) - 1/2, while the value of {dK/dri){l/4) can be 
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FIG. 1: The frequency at the EOB ISCO for binaries having 
spins parallel to i, with mass ratio q = 1x12/ mi and with 
spin-parameter projections onto the direction of L given by 
Xi = X2 = X- As expected, the frequency increases with x 
for a given mass ratio, while for fixed x it increases with q if 
X ^ 0.9, while it decreases with g if x is almost extremal (see 
text for details). 



hopefully determined when more numerical simulations 
with 77 = 0.1-0.25 become available. 

Hereafter, we will use Eq. (|6.1ip and set /<o = 1.4467. 
The latter is determined by requiring that the ISCO fre- 
quency for extreme mass-ratio non-spinning black-hole 
binaries agrees with the exact result of Ref. |4l|, which 
computed the shift of the ISCO frequency due to the con- 
servative part of self-force (see also Ref. [42j where the 
result of Ref. [4l| was compared to the non-spinning EOB 
prediction which resums the function (|5.4ip a la Padc.).^ 



In Fig. [T] we plot the orbital frequency at the ISCO for 
binaries with mass ratio q = 7712 /mi ranging from 10^® 
to 1 and spins aligned with L. In particular, denoting 
by 5i.2 = Xi,2TOi.2 the projections of the spins along the 
direction of i, we consider binaries with Xi = X2 = X- 
We see that the ISCO frequency increases with the mag- 
nitude of the spins x if the mass ratio is fixed, as ex- 
pected from the test-particle case. Also, if the spins 
are kept fixed and small, the ISCO frequency increases 
with the mass-ratio, as it should be to reproduce the 



* We stress that the most general form of K{rf) can include terms 
depending on a? , with a = ISxerrl/^- In particular, a term not 
depending on 77 and proportional to a? could be determined by a 
calculation similar to that in Ref. [4]| , that is by computing the 
shift of the ISCO frequency caused by the conservative part of the 
self-force, for a non-spinning test-particle in a Kerr spacetime. 



FIG. 2: The final spin parameter Xfln as inferred at the EOB 
ISCO, for binaries having spins parallel to L, with mass ratio 
q ~ 7712 /mi and with spin-parameter projections onto the 
direction of L given by Xi = X2 = X- ^s expected, Xfin 
fiattens for large x in the comparable mass case (see text for 
details). 



results of numerical-relativity simulations (see, e.g., the 
non-spinning EOB models of Ref. However, if the 

spins are close to x = 1^ the ISCO frequency decreases 
when the mass ratio increases. This crossover is mir- 
rored by a similar behavior of other quantities evaluated 
at the ISCO — such as the energy, the orbital angular 
momentum, and the coordinate radius — and its physical 
meaning can be explained as follows. When comparable- 
mass almost-extremal black holes merge, the resulting 
black-hole remnant has a spin parameter that is slightly 
smaller than the spin parameters of the parent black 
holes. This is a consequence of the cosmic censorship 
conjecture (see Ref. [46[ and references therein) which 
prevents black holes with spin x > 1 to be formed [13 . 
Therefore, because in the EOB model the position, and 
therefore the frequency, of the ISCO (together with the 
loss of energy and angular momentum during the plunge) 
regulate the final spin of the remnant, and because for 
an isolated black hole the ISCO frequency increases with 
the spin, any EOB model that satisfies the cosmic cen- 
sorship conjecture must have an ISCO frequency that 
slightly decreases with the mass ratio when x ~ 1 • This 
interpretation can be confirmed by computing the final 
spin of the remnant black hole as estimated at the ISCO. 
We have 



Xfin 



•Si + 5*2 -f Lisco 



pbind \2 

-^isco/ 



(6.12) 



which is plotted in Fig. [2] Although the final spin gets 
slightly larger than 1 for high initial spins (because we are 
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FIG. 3: The final spin parameter Xfln as inferred at the EOB 
ISCO for binaries having spins parallel to L, with mass ratio 
q = ni^lmx and with spin-parameter projections onto the di- 
rection of L given by xi = X2 = Xj compared to the remnant's 
final smn parameter predicted by the formula presented in 
Ref. [ijl ("BR09"), which accurately reproduces numerical- 
relativity results. The EOB model and the BR09 formula 
agree when the mass ratio is small (g = 0.1), because the 
emission during the plunge, merger and ringdown is negligible 
in this case. For q — 0.5 and q = 1, there is an offset, because 
the EOB result, at this stage, neglects the gravitational-wave 
emission during the plunge, merger and ringdown (see text 
for details). 



FIG. 4: The mass loss inferred at the EOB ISCO for binaries 
having spins parallel to L, with mass ratio q = -012/ m\ and 
with spin-parameter projections onto the direction of L given 
by Xi ~ X2 = Xi compared to the total mass lost during the 
inspiral, merger and ringdown, as predicted by the formulas 
presented in Ref. gil ("AEI09") and in Ref. ^ ("RIT09"), 
which reproduce numerical-relativity results, although with 
different accuracies because of the different parameter regions 
they cover (see the text for details). The EOB model and the 
AEI09 and RIT09 fits agree when the mass ratio is small 
[q = 0.1), while there is an offset for q = 0.5 and q = 1. 
The reason is that the ringdown emission, which is negligible 
for small mass-ratios, is not taken into account by our EOB 
model at this stage. 



neglecting here the energy and angular momentum emit- 
ted during the plunge, merger and ringdown), the curves 
are remarkably smooth and monotonic (see the corre- 
sponding Fig. 5 of Ref. [2l|) and they flatten at high ini- 
tial spins, as expected. In particular, in Fig.|3]we focus on 
mass ratios g = 1, g = 0.5 and q = 0.1, and plot the final 
spin Xfin as inferred from the ISCO energy and angular 
momentum, together with the final spin of the remnant 
predicted by the formula presented in Ref. [1^ which ac- 
curately reproduce the numerical-relativity results (see 
also Refs. [45l. |49| - |53 | for other formulas for the final spin 
of the remnant). It is remarkable that in spite of the off- 
set between the predictions of the formula of Ref. [i^ and 
the EOB result, which is due to neglecting the energy and 
angular momentum emitted during plunge, merger and 
ringdown, the qualitative behavior of the curves in Fig. [3] 
is the same. Also, we observe that the difference between 
corresponding curves decreases with the mass ratio, with 
the EOB and the numerical-relativity-based results be- 
ing in very good agreement for q = 0.1. This happens 
because the energy and angular momentum emitted dur- 
ing plunge, merger and ringdown become negligible for 



small mass-ratios.^. Similarly, in Fig.|3]we plot the bind- 
ing energy at the ISCO for mass ratios q — 1, q = 0.5 and 
g = 0.1 and compare it with fits to numerical-relativity 
data for the total mass radiated in gravitational waves 
during the inspiral, merger and ringdown. In particular, 
for the q = 1 case we use the fit in Ref. [H ("AEI09"), 
while for q = 0.5 and q = 0.1 we use the fit recently pro- 
posed by Ref. [Hi ("RIT09"). While the AEI fit is more 
accurate than the RIT one for the particular configura- 
tion considered here (see Fig. 11 and related discussion in 
Ref. [3]), the AEI fit is only applicable for comparable- 
mass binaries^^. and for this reason we resort to the more 



® This can be seen by noting that, for a test-particle with mass 
m around a black holes with mass M , the final plunge lasts a 
dynamical time ~ M while the inspiral from large radii to 
the ISCO lasts ~ M^/m. 

The limited applicability of the AEI fit (which is only valid for 
equal-mass binaries with spins aligned or anti-aligned) is indeed 
one reason why it turns out to be more accurate, for the config- 
uration under consideration, than the RIT fit, which is instead 
applicable to more generic binaries. 
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FIG. 5: The frequency at the EOB ISCO for binaries having 
spins parallel to L, with mass ratio q = 1x12/ mi and with 
spin-parameter projections onto the direction of L given by 
Xi = — X2 = X- -A^s expected, the frequency is constant in 
the equal-mass case, because the spins of the two black holes 
cancel out (see text for details). 



FIG. 7; The maximum of the EOB orbital frequency during 
the plunge, for binaries having spins parallel to L, with mass 
ratio q = m2/mi and with spin-parameter projections onto 
the direction of L given by Xi = X2 = X- expected, the 
frequency increases with x for a given mass ratio, while for 
fixed X it increases with q if X ^ 0.9, while it decreases with 
q is X is almost extremal (see text for details). 
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FIG. 6: The final spin parameter Xfln as inferred at the EOB 
ISCO, for binaries having spins parallel to L, with mass ratio 
q = m2/mi and with spin-parameter projections onto the 
direction of L given by xi = — X2 = X- The results are the 
same for all equal-mass binaries, for which the spins of the 
two black holes cancel out (see text for details). 



general RIT fit in the q = 0.5 and q = 0.1 cases. [In Fig.|4] 
we show the predictions of both the AEI and the RIT fit 
in the q = 1 case. Being the AEI fit more accurate, its 
difl:erence from the RIT fit gives an idea of the error bars 
which should be applied to the predictions of the RIT fit 
for q = 0.5 and q = 0.1.] 

In Figs.[5]and[6]we present similar results, for the ISCO 
frequency and for the final spin estimated at the ISCO, 
in the case of spins antialigned with the orbital angular 
momentum. The most apparent feature of these figures 
is that, in the equal-mass case, the quantities under con- 
sideration are independent of x- This happens because 
in this case the spins and S2 are equal and opposite, 
which results in a zero value for the spins S'xcrr and S* 
entering the EOB Hamiltonian (|5.69p . As such, in the 
EOB model, equal-mass binaries with equal and oppo- 
site spins behave as non-spinning binaries. This feature, 
which is also shared by the PN-cxpandcd Hamiltonian, 
until the PN order which is currently known, is also in 
agreement with the results of numerical simulations. In 
fact, equal-mass binaries with equal and opposite spins 
would be indistinguishable with LISA, Virgo and LIGO 
observations [13, H^- Except for this feature, and simi- 
larly to the aligned case discussed above, the behavior of 
the curves in Figs. [S] and [B] is quite smooth and regular 
when going from the equal-mass case to the test-particle 
case. 

In Fig. [7] wc plot the maximum value of the orbital 
frequency during the plunge subsequent to the inspi- 
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FIG. 8: For binaries having spins parallel to L, with mass ratio q — 1x12/ m\ — 1 (left panel) and q — ra2/ra\ =0.1 (right panel) 
and with spin-parameter projections onto the direction of L given by Xi = X2 = Xi ^6 plot twice the maximum of the EOB 
orbital frequency during the plunge against the frequencies of the first 8 overtones of the I = 2, m = 2 quasi-normal mode of a 
Kerr black hole. The quasi-normal mode frequency is computed using the final spin and final mass of the remnant. The final 
spin is estimated by applying the formula of Ref. [i^l, while for the final mass we use the formula of Ref. [iJl ("AEI09") in 
the q = 1 case and that of Ref. [i^ ("RIT09") in the q = 0.1 case. As can be seen, in the g = 0.1 case the peak frequencies 
lie among the high overtones of the i = 2, m = 2 mode, while in the q = 1 case they are generally lower than them. In the 
q — I case we also mark with a square the numerical gravitational-wave frequency at the peak of the /122 mode when x = 0. 
This gravitational-wave frequency coincides with (twice) the maximum of the EOB orbital frequency at the time when the 
matching of the quasi- normal modes is performed in the non-spinning EOB model of Ref. [?[. The numerical gravitational- wave 
frequency is computed from the numerical simulation of Refs. [i^ ("Caltech-Cornell"). 



ral, for binaries with mass ratio q = 7712/7711 and with 
Xi = X2 = X- More precisely, we assume that the particle 
starts off with no radial velocity at the ISCO (thus hav- 
ing angular momentum Lisco and energy i?isco), and 
we compute Pr assuming that the energy and angular 
momentum are conserved during the plunge. Wc find 
that the orbital frequency presents a peak for any value 
of the spins and any mass ratio, and we denote the value 
of the frequency at the peak with Mflj^ax- We note that 
the behavior of M f2niax as a function of the mass ratio is 
similar to that of MfJ^g^Q. In particular, its dependence 
on 77 changes sign when going from small to large spins. 

The physical interpretation of the peak of the orbital 
frequency is that the frequency increases as the effec- 
tive particle spirals in, but when the effective particle 
gets close to the black hole, the orbital frequency has 
to decrease because the particle's motion gets locked to 
the horizon (this is a well-known phenomenon, see for 
instance Ref. [s^, for some of its effect on the test- 
particle dynamics). Said in another way, the orbital fre- 
quency of the effective particle for an observer at infinity 
goes to a constant (or to zero in the non-spinning case Q) 
on the EOB horizon. As a consequence, the peak in the 
frequency can be used to signal the transition between 



two regimes [9[: one in which the deformed black hole 
and the effective particle have different frequencies and 
one in which the two bodies basically move and radiate as 
a single perturbed black hole. For this reason the peak of 
the frequency provides the EOB approach with the nat- 
ural point where to switch to the one-body description, 
i.e., the point where to start describing the gravitational 
waveforms as a superposition of quasi-normal modes. 

We find that the values of MOmax are roughly those 
needed to attach the quasi-normal modes used in EOB 
models to describe the merger and the ringdown ■ 
This is shown in Fig. [51 where we plot twice the maximum 
of the orbital frequency, Mf222 for binaries with q = 1 
(left panel) and q ^ 0.1 (right panel) and with spins 
Xi = X2 = X- We compare MQ22 with the frequency of 
the first 8 overtones of the £ = 2, m = 2 quasi- normal 
mode of a Kerr black hole, computed using the final spin 
and the final mass of the black- hole remnant jsl]. [The 
final spin parameter is estimated by applying the for- 
mula of Ref. |43|. while for the final mass we use the 
formula of Ref. 1441 ("AE109") in the q = 1 case and 
the one of Ref. [U ("R1T09") in the g = 0.1 case.] In 
the 9=1 case we also mark with a square the numerical 
gravitational- wave frequency at the peak of the /122 mode 
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FIG. 9; The maximum of the EOB orbital frequency during 
the plunge, for binaries having spins parallel to L, with mass 
ratio q = milmx and with spin-parameter projections onto 
the direction of L given by xi = — X2 = X- The results are 
the same for all equal-mass binaries, for which the spins of 
the two black holes cancel out (see text for details). 



when X = 0. This gravitational- wave frequency coincides 
with (twice) the maximum of the EOB orbital frequency 
at the time when the matching of the quasi-normal modes 
is performed in the non-spinning EOB model of Ref. . 
The numerical gravitational- wave frequency is computed 
from the numerical simulation of Refs. ("Caltech- 
Cornell"). As can be seen, while in the q = 0.1 case 
the peak frequencies lie among the high overtones of the 
£ = 2, m = 2 quasi-normal mode, in the q — \ case 
they are generally lower than them. Quite interestingly, 
we find that the values of M^^n for x > 0.4 can be in- 
creased up to the frequencies of the quasi-normal modes 
by assuming ijj^2 ^ 30-70 in Eq. (|5.40p . Nevertheless, 
the frequencies that we obtain are comparable to those 
used for the matching with the quasi-normal modes in 
Ref. 0, [13 1 and we therefore expect such a matching to 
be possible also in our EOB model. 

Also, it is interesting to note that the position rmax of 
the frequency peak is quite close (to within 8%) to the 
position of the light ring (or circular photon orbit). This 
fact, which holds exactly in the non-spinning case Q, 
further confirms that the potential barrier for massless 
particles (such as gravitational waves) lies at r ^ J'max- 

Finally, in Fig.[9]wc show the maximum value of the or- 
bital frequency during the plunge for binaries with mass 
ratio q = mijrnx and with xi = — X2 = X- for the 
ISCO quantities, the dependence on the spins and the 
mass ratios is much simpler than in the aligned case, with 
the black-hole spins cancelling out in the equal mass-case 
and thus giving results which arc independent of x- Also, 



we see a smooth transition from the equal-mass to the ex- 
treme mass-ratio case, that our model reproduces exactly. 
Also in this antialigned case, the radius r,„ax agrees with 
the light-ring position to within 4%. 



VII. CONCLUSIONS 

In this paper, building on Ref. [l^l, we computed the 
Hamiltonian of a spinning test particle, at linear order in 
the particle's spin, in an axisymmetric stationary metric 
and in quasi-isotropic coordinates. Then, by applying a 
coordinate transformation, we derived the Hamiltonian 
of a spinning test particle in Kerr spacetime in Boyer- 
Lindquist coordinates 

Wc used those results to construct an improved EOB 
Hamiltonian for spinning black holes. To achieve this 
goal, we followed previous studies [H, [2l| and mapped 
the real two-body dynamics into the dynamics of an ef- 
fective particle with mass /i and spin S* moving in a 
deformed-Kerr spacetime with spin Skcit, the symmet- 
ric mass-ratio of the binary, 77, acting as the deformation 
parameter. 

To derive the improved EOB Hamiltonian, we pro- 
ceeded as follows. First, we applied a suitable canonical 
transformation to the real ADM Hamiltonian and worked 
out the PN-expanded effective Hamiltonian through the 
relation 
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(7.1) 



Then, we found an appropriate deformed-Kerr met- 
ric such that the corresponding Hamiltonian, when ex- 
panded in PN orders, coincided with the PN-expanded 
effective Hamiltonian through 3PN order in the non- 
spinning terms, and 2.5PN order in the spinning terms. 

The (resummed) improved EOB Hamiltonian is then 
found by inverting Eq. (|7.ip . which gives 



improved 
real 




(7.2) 



with Hcs given by Eq. (|5.70p . where a, and 7''' are ob- 
tained by inserting Eqs. (|5.36a|) - (|5.36e|) into Eqs. (|5.44|) - 
(|5.46|) : where Hs is obtained by inserting Eqs. ()5.3ip . 
dSH t, and Eqs. (iCT7a|) - (|E5^ into Eqs. (|i471) . 
and (|4.19p : and where the effective particle's spin S* 
and the deformed-Kerr spin SKcrr (with a = jSKcn lM^) 
arc expressed in terms of the real spins by means of 
Eqs. (it^gl) . (lEMl) . (Ol) and (lOl) . 

The crucial EOB metric potential for quasi-circular 
motion is the potential (r) (which reduces in the non- 
spinning case to the radial potential A{r) of Refs. [9l.[l0j). 
To guarantee the presence of an inner and outer horizons 
in the EOB metric, we proposed to re-write the poten- 
tial At(r) in a suitable way [see Eqs. (|5.7ip and (|5.73p ]. 
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introducing the adjustable EOB parameter Kirj) regu- 
lating the higher-order, unknown PN terms. The rea- 
son why we did not re- write the potential At(r) using 
the Pade summation 21 1 is because Ref. [l^ found that 
when including non-spinning terms at 4PN and 5PN or- 
der, the Pade summation produces spurious poles, does 
not always ensure the presence of an ISCO for binaries 
with spins parallel to L and, even when it docs, the spin 
dependence of physical quantities evaluated at the ISCO 
is quite unusual. 

Restricting the study to circular orbits in the equato- 
rial plane and assuming spins aligned or antialigned with 
the orbital angular momentum, we investigated several 
features of our improved EOB Hamiltonian. Using an 
expression of the EOB adjustable parameter K{ri) which 
reproduces the self-force results in the non-spinning ex- 
treme mass-ratio limit (ill, \^ , we computed the orbital 
frequency at the EOB ISCO, we estimated the final spin 
from the EOB ISCO, and the maximum orbital frequency 
during the plunge. We found that these predictions are 
quite smooth and regular under a variation of rj and of 
the black-hole spins. Quite interestingly, the maximum of 
the orbital frequency during the plunge always exists and 
is close to the light-ring position, as in the non-spinning 
case 0. For this reason, as in the non-spinning case 
the orbital-frequency peak can be used within the EOB 
to mark the matching time at which the merger and ring- 
down start, i.e, the time when, in the EOB formalism, 
the gravitational waveforms start being described by a 
superposition of quasi-normal modes. This will be useful 
in future comparisons of the EOB model with numerical- 
relativity simulations. 

The results of Sec. IVII are an example of the per- 
formances that our improved EOB Hamiltonian can 
achieve. We expect several refinements to be possibly 
needed when comparing our EOB model with accurate 
numerical-relativity simulations of binary black holes. 
We may, for example, extend our model to reproduce also 
the next-to-leading order spin-spin couplings, which are 
known and appear at 3PN order [28l-i30. 33, .36| . Also, we 
might introduce a different mapping between the black- 
hole spins Si, S2 and S*, Sxcrrj a different form of the 
adjustable parameter K(ri), and re- write differently the 
EOB metric potential A((r). We could also introduce 
in it the adjustable parameters 05 and ae at 4PN and 
5PN order, respectively. Moreover, other choices of the 
reference tetrad used to work out the Hamiltonian for 
a spinning test-particle in an axisymmetric stationary 
spacetime could be in principle used, leading to a differ- 
ent (canonically related) EOB Hamiltonian. Lastly, the 
mapping (|7.ip - (|7.2p could me modified by introducing a 
dependence on the spin variables. 

In conclusion, the most remarkable feature of our im- 
proved EOB Hamiltonian is that in the extreme mass- 
ratio limit, it exactly reproduces the Hamiltonian of a 
spinning test particle in a Kerr spacetime, at linear order 
in the particle's spin and at all PN orders. 
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Appendix A: Incorporating spin-spin couplings in 
the effective-one-body Hamiltonian 

The Hamiltonian for a spinning particle in a Kerr 
spacetime that we derived in Sec. IIVI and the Hamilto- 
nian for a spinning particle in a deformed-Kerr spacetime 
that we derived in Sec. IV Cl are only valid at linear order 
in the particle's spin. However, as suggested in Ref. (l9| . 
we can introduce the terms that are quadratic in the par- 
ticle's spin by modifying the quadrupolc moment of the 
Kerr metric. 

In particular, we can add a quadrupolc which is 
quadratic in the particle's spin to the quadrupolc of the 
Kerr metric (which is quadratic in Skctv)- The expres- 
sion for the metric perturbation corresponding to a slight 
change of the Kerr quadrupole can be extracted from 
the Hartle-Thorne metric [s^, H^l, which describes the 
spacetime of a slowly rotating star. Ref. [6l| gives this 
expression in quasi-Boyer-Lindquist coordinates (i.e., in 
coordinates that reduce to Boycr-Lindquist coordinates 
if the quadrupole perturbation is zero, thus reducing the 
spacetime to pure Kerr). This is exactly what is needed 
for our purposes, since we work in quasi-Boyer-Lindquist 
coordinates too. 

In particular, our procedure for introducing the terms 
which are quadratic in the particle's spin into our Hamil- 
tonian consists of modifying the effective metric (|5.36ap - 
(|5.36cp by adding the quadrupolc metric 



1 



j^j4 V -3^ ri , 

where the quadrupole tensor Qij is given by 



^ij — ^ij 37I4 Tlj , 



and h'^'^ is given by [61 1 
1 



1 - 2M/r 
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The functions J^i.2{r) in the above equation are derived 
in Ref. (6l| by transforming the Hartle-Thorne metric to 
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quasi-Boyer-Lindquist coordinates, and are given by 



8Mr{r- 2M) 

15r(r-2M) , / r 
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(2M^ - 3A'/r- Sr^) + 
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Because at large radii Tiir) ~ —Fi[r) « (M/r)^, we see 



that the deformation h'^'^ , when inserted in the Hamil- 
tonian (|5.43p . gives the correct leading-order (2PN) cou- 
pling of the particle's spin with itself. Keeping only the 
leading-order term created by h^'^, the effective EOB 
Hamiltonian therefore becomes 
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Hs + /3'p^ + a ^m:^ p,pj + Qi{p) 
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